*use corruption_scandal.dta

keep aware lrself like_party1 news_watch5 dif_news generation female education ob_incom HN YN bi_bh_aide_great_fault only_park_fault bi_opp_responsible no_pro_impeach2 no_anti_impeach2 near_seoul c_hobby_meeting c_help_meeting bi_no_anti_impeach2 ID



*** Table 3 Results 


oprobit aware lrself like_party1 news_watch5 dif_news, nolog  

oprobit aware lrself like_party1 c.news_watch5##i.dif_news, nolog  

oprobit aware lrself like_party1 news_watch5 dif_news generation female education ob_incom HN YN, nolog 

oprobit aware lrself like_party1 c.news_watch5##i.dif_news generation female education ob_income HN YN, nolog   


margins if dif_news==1, at (news_watch5=(1(1)5)) vsquish post predict(outcome(4))
marginsplot, recast (line) recastci(rarea) name(a, replace)

oprobit aware lrself like_party1 c.news_watch5##i.dif_news generation female education ob_income HN YN, nolog   


margins if dif_news==0, at (news_watch5=(1(1)5)) vsquish post predict(outcome(4))
marginsplot, recast (line) recastci(rarea) name(b, replace)


graph combine a b 


** Table 4 Results 

probit bi_bh_aide_great_fault news_watch5 dif_news like_party1 lrself generation female education YN HN ob_income, nolog  
 
probit bi_bh_aide_great_fault c.news_watch5##i.dif_news like_party1 lrself generation female education YN HN ob_income, nolog 

probit only_park_fault news_watch5 dif_news like_party1 lrself generation female education YN HN ob_income, nolog 

probit only_park_fault c.news_watch5##i.dif_news like_party1 lrself generation female education YN HN ob_income, nolog 

probit bi_opp_responsible news_watch5 dif_news like_party1 lrself generation female education YN HN ob_income, nolog

probit bi_opp_responsible c.news_watch5##i.dif_news like_party1 lrself generation female education YN HN ob_income, nolog  
 

 ** Table 5 Results 
 

nbreg no_pro_impeach2 news_watch5 dif_news lrself like_party1 generation female education HN YN near_seoul ob_income, nolog 


nbreg no_pro_impeach2 c.news_watch5##i.dif_news lrself like_party1 generation female education HN YN near_seoul ob_income, nolog 



zinb no_anti_impeach2 news_watch5 dif_news lrself like_party1 generation female education HN YN near_seoul ob_income, nolog inflate (ID) zip


zinb no_anti_impeach2 c.news_watch5##i.dif_news like_party1 lrself  generation female education YN HN near_seoul ob_income, nolog inflate (ID) zip 


margins if dif_news==1, at (news_watch5=(1(1)5)) vsquish post 
marginsplot, recast (line) recastci(rarea) name(c, replace)


zinb no_anti_impeach2 c.news_watch5##i.dif_news like_party1 lrself  generation female education YN HN near_seoul ob_income, inflate (ID) zip 


margins if dif_news==0, at (news_watch5=(1(1)5)) vsquish post 
marginsplot, recast (line) recastci(rarea) name(d, replace)

graph combine c d


** Table 6 Results 

cmp(aware=news_watch5 like_party1 dif_news generation female education YN HN lrself ob_income)(news_watch5=aware like_party1 dif_news generation female education c_hobby_meeting c_help_meeting lrself YN HN ob_income), nolog nolr ind($cmp_oprobit $cmp_oprobit)

cmp(bi_bh_aide_great_fault=news_watch5 dif_news like_party1 lrself generation female education YN HN ob_income)(news_watch5= bi_bh_aide_great_fault like_party1 generation female education dif_news c_hobby_meeting c_help_meeting YN HN lrself ob_income), nolog nolr ind($cmp_probit $cmp_oprobit)

cmp(only_park_fault=news_watch5 dif_news like_party1 lrself generation female education YN HN ob_income)(news_watch5= only_park_fault like_party1 generation female education dif_news c_hobby_meeting c_help_meeting YN HN lrself ob_income), nolog nolr ind($cmp_probit $cmp_oprobit)

cmp(bi_opp_responsible=news_watch5 dif_news like_party1 lrself generation female education YN HN ob_income)(news_watch5= bi_opp_responsible like_party1 generation female education dif_news c_hobby_meeting c_help_meeting YN HN lrself ob_income), nolog nolr ind($cmp_probit $cmp_oprobit)

cmp(no_pro_impeach2=news_watch5 dif_news like_party1 lrself generation female education near_seoul ob_income)(news_watch5= no_pro_impeach2 like_party1 generation female education dif_news c_hobby_meeting near_seoul lrself ob_income), nolog nolr ind($cmp_cont $cmp_oprobit)
 
cmp(bi_no_anti_impeach2=news_watch5 dif_news like_party1 lrself generation female education HN YN near_seoul ob_income)(news_watch5= bi_no_anti_impeach2 like_party1 generation female education dif_news c_hobby_meeting c_help_meeting HN YN near_seoul lrself ob_income), nolog nolr ind($cmp_probit $cmp_oprobit)





